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More on SU(3) lattice gauge theory in the fundamental-adjoint plane* 
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We present further evidence for the bulk nature of the phase transition line in the fundamental-adjoint action 
plane of SU(3) lattice gauge theory. Computing the string tension and some glueball masses along the thermal 
phase transition line of finite temperature systems with Nt = 4, which was found to join onto the bulk transition 
line at its endpoint, we find that the ratio y/o '/T c remains approximately constant. However, the mass of the 
0++ glueball decreases as the endpoint of the bulk transition line is approached, while the other glueball masses 
appear unchanged. This is consistent with the notion that the bulk transition line ends in a critical endpoint with 
the continuum limit there being a 4> l theory with a diverging correlation length only in the ++ channel. 



1. INTRODUCTION 

Studies in the early 80 's of the phase diagram 
of fundamental-adjoint pure gauge systems, 

S = Pi D 1 - jjKeTrUp] 
p 

+ p°T,t L -jp TiU P TiUp i' w 
p 

revealed a non-trivial phase structure with first 
order (bulk) transitions in the region of small f3j 
[1,2]. In particular a line of first order transitions 
points towards the fundamental axis and, after 
terminating, extends as a roughly straight line 
of bulk crossovers to the fundamental axis and 
beyond. 

In a couple of recent papers Gavai, Grady and 
Mathur [3] and Mathur and Gavai [4] returned to 
investigating the behavior of pure gauge SU(2) 
theory in the fundamental-adjoint plane at fi- 
nite, non-zero temperature. They raised doubts 
about the bulk nature of the phase transition and 
claimed that their results were consistent with 
the transitions, for lattices with temporal extent 
Nt = 4, 6 and 8, to be of thermal, deconfining 
nature, displaced toward weak coupling with in- 
creasing Nt ■ On the transition line for a fixed Nt 
there is then a switch from second order behavior 
near the fundamental axis to first order behavior 
at larger adjoint coupling. In a Landau Ginzburg 
model of the effective action in terms of Polyakov 
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lines, Mathur [5] reported that he could reproduce 
the claimed behavior seen in the numerical simu- 
lations. These results, should they be confirmed, 
are rather unsettling, since they contradict the 
usual universality picture of lattice gauge theory 
with a second order deconfinement transition for 
gauge group SU(2). 

Puzzled by the finding of Ref. [3], we stud- 
ied the finite temperature behavior of pure gauge 
SU(3) theory in the fundamental-adjoint plane 
[6,7]. We obtained results in agreement with the 
usual universality picture: there is a first order 
bulk transition line ending at 

(P* S ,K) = (4.00(7), 2.06(8)). (2) 

The thermal deconfinement transition lines for 
fixed Nt (being of first order for gauge group 
SU(3)) in the fundamental adjoint plane are or- 
dered such that the thermal transition for smaller 
Nt occurs to the left, at smaller f3j, than that 
for a larger Nt- In this order the thermal transi- 
tion lines join on to the bulk transition line. The 
thermal transition line for Nt = 4 joins the bulk 
transition line very close to the critical endpoint. 
This is shown in Figure 4 of Ref. [7], reproduced 
here as Figure 1. 

To solidify this picture, which is in agree- 
ment with the usual universality scenario, we 
have continued the investigation studying zero- 
temperature observables. We computed the 
string tension and the masses of some glueballs, 
in particular the ++ glueball, along the thermal 
transition line for Nt = 4. For universal contin- 
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Figure 1. The phase diagram together with the 
thermal deconfinement transition points for Nt = 
2, 4, 6 and 8 from Ref. [7]. The lower plot shows 
an enlargement of the region around the end point 
of the bulk transition. 

uum behavior and the glueball masses should 
be constant along the thermal transition line for 
a fixed Nf, leading to constant ratios T c /~Ja and 

Of course, at small N t , corresponding to a large 
lattice spacing a, we expect to see some devia- 
tions from this constant behavior. However, we 
find large deviations for m ++ / as we approach 
the critical endpoint of the bulk transition along 
the Nt = 4 thermal transition line. The scalar 
glueball mass decreases dramatically; much more 
than what could be expected from simple scaling 
violations at a large lattice spacing. On the other 
hand, this is not really a surprise since at the crit- 
ical endpoint at least one mass in the ++ channel 
has to vanish. We elaborate on our findings in the 
next sections and then discuss the implications 
for the scaling behavior along the fundamental 
(or Wilson) axis. 

2. OBSERVABLES AND ANALYSIS 

We have made simulations of the model with 
action (1) along the thermal transition line for 



Nt = 4, and continued along the bulk transition 
line, on a 12 4 lattice. We computed finite T ap- 
proximants to the potential from time-like Wil- 
son loops constructed with 'APE'-smeared spa- 
tial links [8] to increase the overlap to the ground 
state potential. On and off axis spatial paths were 
considered with distances R = n, \/2n, \/%n and 
\/bn, with n = 1, 2, ...an integer. The string 
tension was then extracted from the usual fit 

V(R) = Vo~^ + l (G L (R) - £) + aR. (3) 

Here Gl(R) is the lattice Coulomb potential, in- 
cluded in the fit to take account of short distance 
lattice artefacts. Our fits are fully correlated \ 2 ~ 
fits with the correlations estimated by bootstrap, 
after binning to alleviate autocorrelation effects. 
The results of the best fits are listed in Table 1. 

We also computed glueball correlation func- 
tions in the ++ , 2 ++ and l" 1 channel that 
can be built from simple plaquette operators, but 
built from the smeared links, already used for the 
computation of the potential. Not surprisingly 
for computations done around the critical cou- 
pling for the N t = A thermal phase transition, 
we did not obtain a significant signal in the l" 1 
channel and only an effective mass from time dis- 
tances t = 0/1 in the 2 ++ channel - we had 500 
measurements everywhere, except for f3 a = 2.0, 
near the critical endpoint where the number was 
increased, as given in Table 1. In the ++ channel 
we got a signal at distance t = 1/2 at small f3 a 
and out to t = 3/4 at f3 a = 2.0. Our best results 
are also given in Table 1. 

The quantities ^/a and m ++ are shown in Fig- 
ure 2 plotted versus f3 a . As can be seen, ^/a re- 
mains approximately constant along the thermal 
transition line for Nt = 4 - the errors shown in 
the figure are statistical only; no estimate of the 
error from the uncertainty in the determination 
of /3f c has been attempted except for f3 a = 2.0. 
There, the computation has been repeated for two 
nearby couplings, also listed in Table 1; the vari- 
ation with /3f becomes so rapid that the error in 
the determination of f3j c becomes the dominating 
factor. 

While our estimate for m 2 ++ , albeit an unreli- 
able estimate since we had to use distances t = 
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Pa 


Pjc(N t = 4) 


Pj 


L 


N meas 




m ++ 


m 2 ++(t = 1) 


0.0 


5.6925(2) 


5.7 a 






0.4099( 12) 


0.97( 2) 


2.39(13) 


0.5 


5.25(5) 


5.25 


12 


500 


0.4218( 28) 


0.93(11) 


2.29(13) 


1.0 


4.85(5) 


4.85 


12 


500 


0.4024( 82) 


0.78(28) 


2.45(19) 


1.5 


4.45(5) 


4.45 


12 


500 


0.3743( 51) 


0.56(17) 


2.13( 9) 


2.0 


4.035(5) 


4.03 


12 


1000 


0.555 ( 11) 


0.34( 6) 


3.11(16) 


2.0 


4.035(5) 


4.035 


12 


2000 


0.4725(128) 


0.20( 4) 


3.17(14) 


2.0 


4.035(5) 


4.04 


12 


1000 


0.3750( 24) 


0.27( 8) 


2.37( 9) 


2.25 


3.8475(25) 


3.8475' 


12 


500 


0.619 ( 18) 


0.37(10) 


3.15(37) 


2.25 


3.8475(25) 


3.8475 c 


12 


500 


0.3005( 22) 


0.61( 8) 


1.59( 7) 


2.25 


3.8475(25) 


3.8475 c 


16 


500 


0.2965( 19) 


0.62( 5) 


1.72( 9) 



Table 1 

The results in the neighborhood of the Nf = 4 thermal transition line. Comments: (a) at f3 a = 0.0 
comes from [9], m ++ from [10]; we did not take their best value, but rather the effective mass from the 
same distance as at f3 a = 0.5; m 2 ++ is from [11]; (b) in the disordered phase and (c) in the ordered phase 
on the bulk transition. 



and 1 to obtain enough of a signal to extract an 
effective mass, also remains approximately con- 
stant, m ++ shows a remarkable decrease as the 
critical endpoint is approached. This observed 
behavior suggests that the mass in the ++ chan- 
nel vanishes at the critical endpoint, thereby giv- 
ing strong additional evidence for the existence of 
this critical endpoint, since at a critical point at 
least one mass, in the ++ channel to have a ra- 
tionally invariant continuum limit, must vanish. 

Since no other observable seems to be dramati- 
cally affected by the critical endpoint, we conjec- 
ture that the continuum theory one would obtain 
there is simply the (trivial) (j> 4 theory. To sub- 
stantiate this claim somewhat, we made a fit to 
the scalar mass of the form 

m ++ = A(p* a -p a ) p (4) 

expected near a critical point. A 3-parameter fit 
gave A = 0.76(11), p = 0.35(20) and /?* = 2.02(6) 
with a x 2 of 0.29 for 2 dof. Note that the estimate 
for /?* is in agreement with the previous estimate 
(2) obtained in [7] from fits to the jump in the 
plaquette across the bulk transition line. Within 
its large error, the exponent p is compatible with 
the mean field value 0.5 of (j> 4 theory, up to log- 
arithmic corrections. Since the errors of the fit 
parameters are rather large we also made a fit 
with /?* held fixed at its value 2.06 of (2). This 



fit gave A = 0.71(3) and p = 0.44(5) with a \ 2 
of 0.47 for 3 dof. Again, p is compatible with the 
mean field value. Finally, a fit with /?* = 2.06 
and p = 0.5 both held fixed gave A = 0.68(1) 
with a x 2 °f 1-50 for 4 dof. This last, still very 
acceptable fit is included in Figure 2. 

3. IMPLICATIONS FOR SCALING 

In the previous section we have corroborated 
the existence of a first order bulk phase transi- 
tion line ending in a critical endpoint. We have 
provided evidence that physical observables are 
little affected by this phase transition line. The 
notable exception to this is the glueball mass in 
the ++ channel which decreases as the critical 
endpoint is approached and vanishes there. 

The influence of the critical endpoint on the 
++ glueball appears still visible in the crossover 
region on the fundamental (Wilson) axis. This 
has first been argued in Ref. [12]. The ++ glue- 
ball is lighter in the crossover region, leading to a 
different scaling behavior than other observables. 
This can be seen in Figure 3 where we show the 
latest data of T c /^/a from Ref. [13] and m ++ / 
with ^/a taken from Refs. [9,14] and the glueball 
mass from Refs. [10,15,16]. While T c /^/a stays 
approximately constant in the [3 interval shown, 
m ++ 1 decreases visibly in the crossover re- 
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Figure 2. ^/a (octagons) and m ++ (squares) 
as a function of f3 a along the thermal transition 
line for N t = 4. At (3 a = 2.25 we show the 
results from both phases at the bulk transition. 
The dashed vertical line gives the location of the 
critical endpoint, (2), with the dotted lines in- 
dicating the error band. The curve is a fit to 
m 0++ = A{2M- pa) 1 ' 2 . 

gion around [3 = 5.7. 

It has long been known that the scalar glue- 
ball mass scales differently in the crossover re- 
gion than T c and the string tension. The scalar 
glueball seemed more compatible with asymp- 
totic scaling. However, in view of our findings 
that the different behavior of the scalar glueball 
mass comes from the influence of the nearby criti- 
cal endpoint this asymptotic scaling behavior ap- 
pears to be accidental. 
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